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ABSTRACT 



We present the analysis result of a gravitational binary-lensing event OGLE-2005- 
BLG-018. The light curve of the event is characterized by 2 adjacent strong features 
and a single weak feature separated from the strong features. The light curve exhibits 
noticeable deviations from the best-fit model based on standard binary parameters. 
To explain the deviation, we test models including various higher-order effects of the 
motions of the observer, source, and lens. Prom this, we find that it is necessary to 
account for the orbital motion of the lens in describing the light curve. From modeling 
of the light curve considering the parallax effect and Keplerian orbital motion, we are 
able to measure not only the physical parameters but also a complete orbital solution 
of the lens system. It is found that the event was produced by a binary lens located in 
the Galactic bulge with a distance 6.7 it 0.3 kpc from the Earth. The individual lens 
components with masses 0.9 ±0.3 Mq and 0.5 ±0.1 Mq are separated with a semi-major 
axis of a = 2.5 ± 1.0 AU and orbiting each other with a period P = 3.1 ± 1.3 yr. The 
event demonstrates that it is possible to extract detailed information about binary lens 
systems from well-resolved lensing light curves. 

Subject headings: binaries: general - gravitational lensing: micro 



1. Introduction 

Microlensing occurs when an astronomical object approaches close to the line of sight toward a 
background star. Due to the gravity of the intervening object (lens), light rays from the background 
star (source) bend, causing splits and distortions of the source star image. For Galactic microlensing 
events, the separation between the split images is an order of milli-arcsec and thus it is difficult to 
directly observe the split images. However, the lensing phenomenon can be observed by the change 
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of the source star brightne ss. For a point-s ource single-lens event, the magnification of the source 
star flux is represented by (|Paczvhskilll986l ) 
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(1) 



where u is the lens-source separation in unit of the angular Einstein radius 9e, tp is the time of 
closest lens-source approach, uq is the lens-source separation at that moment, and tE is the time 
required for the source to transit 6e (Einstein time scale). The Einstein radius is related to the 
physical parameters of the lens system by 



Oe = (KM7^rel)^/^ 



(2) 



where k = 4G/ (c^AU), M is the mass of the lens, tTi-cI = A\J{D^^ — Dg^), and Dl and Ds are the 
distances to the lens and source, respectively. A standard single-lens event is characterized by a 
non-repeating, smooth, and symmetri c light curve and mo delin g it requires 3 parameters of to, uq, 
and tE- Since the first detections bv lAlcock et al.l (|l993l ) and lUdalski et al.l (jl993l ) , microlensing 



events have beeri dete cted toward various star fields includi ng the Galactic bulge (lUdalski et al 



2005 
M31 d 



Sumi et al 



2010ll. Large and Small Magellanic Clouds (IWvrzvkowski et al.ll2009l . l20ld ). and 



Calchi Novati et al.ll2010l ). Currently, events are being detected with a rate of more than 500 



events per year, mostly toward the bulge field. 

Among lensing events, a fraction of events are produced by lenses composed of two masses. 
These binary-lens events can exhibit ligh t curves that are dramatically different from those of 



single-lens events (jMao &: Paczvhskilll99ll ). The most prominent features occur when the source 
closely approaches or crosses caustics, which represent the set of source positions at which the 
lensing magnification of a point source becomes infinity. Describing a standard binary-lens light 
curve requires to include three additional parameters: the mass ratio of the companion to its host, 
q, the projected separation between the lens components in units of the Einstein radius, s±, and 
the angle between the source trajectory and the binary axis, a. For a caustic-crossing event, an 
additional parameter of the source radius normalized by the Einstein radius , p^, (normalized source 
radius), is required to account for the finite-sou rce effect (jPominkl Il995l : iGaudi &: Gouldl Il999l : 
Gaudi k Petterslbood : [Peicha Hevrovskvl bood ) . 



However, modeling binary-lens light curves with the standard parameters is occasionally not 
adequate to precisely describe light curves. This is because light curves are subject to various 
higher-order effects that result in deviations from the canonical form. The causes of these effects 
include the motions of the observer, source, and lens during the event. The change of the observer's 
position induced by the orbital motion of the Earth around the Sun causes the source motion with 
respect to the lens to deviate from rectilinear, and thus the resulting light curve ca n exhibit long- 
term deviations. This is known as the "parallax" effect ([Gouldl Il992l : lRefsdallll966l ). If the source 
star is a binary, the source trajectory can also be affected by the orb ital motion of the s ource 
over the course of the event. This is known as the "xallarap" effect (jHan &: Gouldl 119971 ). i.e.. 
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Fig. 1. — Light curve of the microlensing event OGLE-2005-BLG-018. Also presented our light 
curves from modelings based on the standard binary parameters (dotted curve) and considering 
the parallax and orbital effects (solid curve). 



reverse order of the spell of "parallax". Finally, the binary nature of the lens implies that the 
positions of the lens components vary in time due to their orbital motion. This "orbital" effect 
causes not only the change of the source position with respect to th e lens but also t he magnification 
pattern because the projected binary sepa, ration changes in time (jDominikl Il998l : lloka et al.l 11993 : 
Albrow et al.ll200d : iRattenburv et al.ll2002l ). The deviations in lensing light curves caused by these 
second-order effects are usually very small and thus difficult to measure. However, when they are 
measured, they provide information that allows to better constrain the lens system. For example, 
if the deviation caused by the parallax effect is measured, it is possible to determine the p hysical 
param eters of the lens mass, distance to it, and the projected separation in physical units () Gould 
19921 ). If the orbital effect is measured, one can further constrain the lens syst em by determining 



the orbital parameters and intrinsic s eparation between the lens components (iGaudi et al 



Dong et al. 



200S 



Bennett et al.ll2010l : iPennv et al 



201C 



Yee et al 



20101; 



Skowron et al 



2008 



201 



3). 



In this paper, we analyze the light curve of the binary-lensing event OGLE-2005-BLG-018 
by combining all available data obtaine d from survey and follow-up observations. The event was 
analyzed before by lSkowron et al.l (120071 ) based on the OGLE data and their model shows noticeable 
deviations. Our preliminary modeling also indicated that a model based on the standard binary 
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parameters is not adequate enough to precisely describe the observed hght curve, suggesting the 
need to consider higher-order effects. We conduct analysis of the light curve considering various 
second-order effects and present the constraints of the lens system obtained from modeling. 



2. Observation 



The event OGLE-2005-BLG-018 occurred on a Galactic bulge star located at {a, (5)j2000 = 
(17'^51™23'^.53, -29°39'22".8), which corresponds to the Galactic coordinates ( l b) = (0.03°, -1.45°) 



The event was detected by the Optical Gravitational Lensing Experiment (jOGLE: Udalski et al 



20031 ) using the 1.3 m Warsaw telescope of Las Campanas Observatory in Chile. An anomaly 
alert was issued on 2005 Mar 31 by the OGLE group. In addition, a series of real-time models 
were issued by M. Dominik (2010, private communication). Following the alert and models, the 
event was intensively observed by foll ow-up grou ps including the Pro bing Lensing Anomalies Net- 
work dPLANET: Beaulieu et aP bood). RoboNet f|Tsapras et ahlbood ). and Microlensing Follow-Up 
Network (j^FUN: Gouldll2006l ) by using six telescopes located on four different continents. The tele- 
scopes used for follow-up observations include 2.0 m Faulkes Telescope N. (FTN) in Hawaii, 2.0 m 
Liverpool Telescope (LT) in La Palma, Spain, 1.0 m of Mt. Canopus Observatory in Australia, and 
1.54 m Danish Telescope of La Silla Observatory in Chile, 0.6 m of Perth Observatory in Australia, 
and 1.3 m SMARTS telescope of CTIO in Chile. Thanks to the follow-up observations, the light 
curve was very densely resolved. 

Figure [1] shows the light curve of the event. It is characterized by three distinctive features, 
occurring at HJD ~ 2453460 (h), 2453512 (^2), and 2453528 (ts). The two adjacent peaks at ^2 and 
ts are strong while the other peak at ti is relatively weak and separated from the strong features. 



3. Modeling of Second-Order Effects 

We first conduct modeling of the light curve with the set of standard binary parameters. As 
shown in Figure [H the best-fit light curve from this initial modeling shows noticeable deviations 
from the observed light curve especially near the part of the light curve around the weak feature 
although the model light curve describes the two strong features relatively well. Investigation of the 
lens system geometry obtained from the standard modeling indicates that the projected separation 
between the binary lens components is smaller than the Einstein radius, i.e., s_l < 1. In this 
case, the resulting caustics are composed of 3 segments, where one large central caustics is located 
around the cent er of mass of the binary an d the other two small caustics are located away from the 
central caustic ( Scheneider &: Weiss 1986l )FI The model also indicates that the two strong features 



^In more precise term, the number of caustic is 3 when sx < (1 + q')^^*(l + q^^^) and the three caustics merge 
into a single one as the separation becomes equivalent to the Einstein radius. 
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of the light curve were produced by two successive crossings of the source over the central caustic 
and the weak feature was produced by the approach to one of the two peripheral caustics. Events 
produced by such a lens system are susceptible to the orbital effect because the peripheral caustic 
moves considerably even for a small shift of the binary axis. In addition, the long duration of the 
event, which lasted ~ 200 days, raises the need to consider both the parallax and xallarap effects. 
We, therefore, conduct modelings of the light curve considering these second-order effects as well. 

Considering the parallax effect into modeling requires to include two additional parameters of 
tte^n and tte.e- These parameters represent the two components of the lens parallax vector tte 
projected on the sky in the north and east equatorial coordinates, respectively. The direction of 
the parallax vector corresponds to that of the lens-source relative motion in the frame of the Earth 
at a specific time of the event. We set the reference time at the moment of the second perturbation 
peak, i.e., t2- The size of the parallax vector corresponds to the ratio of the Earth's orbit to the 
Einstein radius projected on the observer observer's plane, i.e., tte = A\J /[r^Ds/{Da, — Dl)]- 

Under the assumption of a circular orbit and a very faint binary companion, the xallarap effect 
is described by 5 parameters. They are the orbital period of the source, Ps, inclination, is, phase 
angle, ip, and the two components of the xallarap vector in the north and east direction, S,e,n and 
Ce,e- The magnitude of the xallarap vector corresponds to the ratio of the source star's orbit 
to the Einstein radius projected on the source plane. 

To account for the orbital effect, we consider 2 types of parameterization. The first one is based 
on the approximation that the binary lens rotates with a constant angular speed and the projected 
separation between its components changes with a constant rate. This parameterization requires 
two parameters of da/dt and ds±/dt, which represent the changing rates of the angle between the 
binary axis and the source trajectory and the projected separation between the lens components, 
respectively. In the second parameterization, we fully consider the Keplerian orbital motion. This 
requires to include 2 extra parameters in addition to the orbital parameters used in the first type 
of parameterization. These additional parameters are sy and ds\\/dt, where sy represents the line- 
of-sight separation between the binary components in units of 0e and ds\\/dt represents its rate of 
ch ange. For the full des cription of the orbital lensing parameters, see the summary in the Appendix 



of ISkowron et all mm . 



With these parameterizations, we test 6 different models. The first model is based on standard 
binary parameters with a static lens and source and no parallax motion of the Earth (standard 
model). The second and third models include the parallax and xallarap effects, respectively. The 
fourth model includes the orbital effect of the lens with 2 parameters of da/dt and ds±/dt. In 
the last two models, we consider both the parallax and orbital effects where the orbital effect is 
described by 2 and 4 parameters, respectively. We note that 4-parameter orbital modeling must 
include the parallax parameters whereas 2-parameter orbital modeling does not necessarily needs 
this. 



We search for the solution of the best-fit parameters of the individual models by minimizing 
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Fig. 2. — Residual of data from 6 different models. 



in the parameter space. For binary-lensing modeling, this is a very complex procedure due to several 
reasons. First, the complexity of the surface i n the paramete r space makes it difficult to rule out 
the possibility of the existence of local minima (jDominiklll999l ). implying that even if a plausible 
model is found, it is difficult to be sure the solution is the correct one. This makes it difficult to 
use a simple downhill approach to search for solutions. Second, modeling is further complicated 
by the sheer size of the parameter space. The large number of parameters implies that brute-force 
searches for solutions are very difficult and extremely time-consuming. To resolve the degeneracy 
problem but avoiding searches throughout all parameter space, we use a hybrid approach in which 
grid searches are conducted over the space of a set of parameters and the remaining parameters 
are searched by using a downhill approach. We choose s±, q, and a as the grid parameters because 
they are related to the light curve features in a complex way such that a small change in the values 
of the parameters can lead to dramatic changes in the resulting light curve. On the other hand. 
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the other parameters are more directly related to the light curve features and thus they can be 
searched for by using a downhill approach. For the minimization in the downhill approach, we 
use the Markov Chain Monte Carlo (MCMC) method. 

Another difficulty in binary-lensing modeling arises due to large computation. Most binary-lens 
events exhibit perturbations induced by caustic crossings or approaches during which the finite- 
source effect is important. Calculating finite magnifications involves a numerical method, which 
requires heavy computations. Considering that modeling requires to produce a large number of light 
curves of trial models, it is important to apply an efficient method for magnification calculations. 
We accelerate the computation first by minimizing computations in the numerical method and 
second by restricting the numerical computation only when it is necessary. The numerical method 
applied for finite magnification computations is based on the ray-shooting method. In this method, 
a large number of rays are uniformly shot from the image plane, bent according to the lens equation, 
land on the source plane, and then the magnification is computed as the ratio of the number density 
of rays on the source surface to the density on the image plane. In this process, we reduce the 
number of rays required for magnification computations by shooting only ones arriving at the 
region around the caustics. We further restrict numerical co mputations by applying a simple 



analy tic hexadecapole approximation for finite magnifications (|Peicha &: Hevrovskv 
2008 ) unless the source is located very close to the caustics. 



200S 



Gould 



We incorporate the effect of the limb-darkening of the source star surface when we compute 
the finite-source magnification. The surface brightness is modeled by 



1 



1 — - COS ( 



(3) 



where Tx is the linear limb-darkening coefficients. Fx is the flux, and 9 is the angle between the 
normal to the source star's surface and the line of sight toward the star. From the color of the 
source star measured from the location on the color-magnitude diagram, it is found that the source 
is a clump giant. We, therefore, use the limb-darkeni ng coefficients of Fy = 0.708, Tj = 0.613, and 
F/j = 0.508 by adopting the values from lClaretl (|2000l ) under the assumption that i^turb = 2 km s^^, 
log(5/50) = -1.9, and T,s = 4750 K. 



4. Result 

In Table [H we present the solutions found for the 6 tested models. In Figure El we also present 
the residuals of the data from the best-fit light curves of the individual models. We find that 
neither the parallax nor the xallarap effect alone is enough to precisely describe the observed light 
curve although the models considering these effects improve the fit from the standard model with 
Ax^/dof ~ 1113/1600 and 1151/1597, respectively. For both models, the fits are still poor near 
the weak feature of the light curve. 

On the other hand, it is found that the light curve is well described by the model including 
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Fig. 3. — Scatter plot of Markov chains in the space of the parallax parameters tte^b and vte^at for 
the best-fit model. Different colors of points represent Ax^ from the minimum. 



the orbital effect. We find that the fit with 2 orbital parameters is better than the standard model 
by Ax^/dof ~ 1853/1600. It is also better than the best-fit parallax and xallarap models by 
Ax^/dof ~ 740/1600 and 702, respectively. Therefore, we find that the dominant second-order 
effect for the deviation between the observed data and standard model is the orbital motion of the 
lens. As mentioned in the previous section, the importance of the orbital effect was expected due 
to the specific geometry of the lens system in which the source trajectory passes the central and the 
peripheral caustics of a close binary. In this sense, the event has much in common with the event 



MACHO 97-BLG-41 for which the orbital effect was measured for the first time (Albrow et al 



2000) 



The Einstein radius is measured from the normalized source radius, p^,, determined from mod- 
eling combined with the angular radius of the source star, 9^, i.e., 9e = Oi^/p^,. We measure the 
angular source radius by first measuring the de-reddened color of the so urce star and convert ing it 
into radius by using the relation between the color and angular radius of lKervella et al.l (120041 ) . For 
the calibration of the magnitude and color of the source star, we use the centroid of bulge clump 
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Fig. 4. — Histograms of lens orbital parameters obtained from the 4-parameter orbital fit. 



giants in the color-magnitude diagram as a reference u nder the assump tion that the source and 
clump giants experience the same amount of extinction (|Yoo et al.ll2004l ). The measured Einstein 
radius is Be = 0.506 ± 0.044 mas. 

Although the importance of the orbital motion of the lens in describing the observed light 
curve is identified, we conduct additional modeling considering both orbital and parallax effects in 
order to see the possibility of further improvement of the fit and to constrain physical parameters 
of the lens system. The orbital effect is considered by both models with 2 and 4 parameters. Prom 
the model with the parallax plus 2 orbital parameters, it is found that adding the parallax effect 
does not improve the fit significantly. This could be because the parallax is poorly constrained 
or because it constrained close to zero. From the other modeling with 4 orbital parameters, it is 
found that the latter holds. In Figure [3] and HI we present the scatter plot of Markov chains in the 
'^E,E-'^E,N parameter space and the histograms of the microlens orbital parameters, respectively. 
They show that the parallax and orbital parameters are reasonably well constrained. We measure 
the lens parallax of tte = 0.028 it 0.010. A small parallax value suggests that either the lens is heavy 
or it is located away from the Earth. On the top of the light curve in Figure [H we present the 
best-fit light curve for this model. In Figure [5l we also present the source trajectory with respect 
to the caustic. We note that the caustic shape varies with time. We present three sets of caustics 
corresponding to the times of ti, t2, and t^. Also marked are the positions of the lens components 
at the corresponding times. 
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We determine the physical and orbital parameters of the lens system based on the measured 
lensing parameters. This requires to adopt a value of the Einstein radius in the modeling process 
including 4 orbital parameters. We adopt this value as the one measured from the model with 
2 orbital parameters, i.e., 9e = 0.506 ± 0.044 mas. In principle, the value of 6e could change as 
more parameters are added. However, this change is usually very small because the constraint of 
p^, from which 9^ is measured, is provided by the very localized region of the light curve where 
the finite-source effect is important, while the orbital and parallax effects are constrained from 
the overall shape of the light curve. As a result, the physical and orbital parameters are barely 
affected by the adopted value of 9e- With the Einstein radius and the lens parallax determined 
from modeling, the mass and distance to the lens are determined by 

9v 

M = Mi + M2 = — (4) 

KITE 

'''''^ AU 

^ 7Te9e + vts ' ^ ^ 

respectively, where Mi and M2 are the masses of the heavy and light components of the lens, 
respectively, and vrs = AU/Ds is the parallax of the source star. 

In addition to the constraints provided by the light curve itself, the lens system can also be 
constrained by the blended flux. This is because the flux from the lens cannot exceed the measured 
blended flux. We find that the blended flux is negligible compared to the flux from the source star. 
Even considering that the source is a giant, this provides the constraint that the primary of the 
lens should be a main-sequence star. Therefore, we set the upper mass limit of the primary as 
~ 1.3 Mq, and thus the total mass of the lens should be < 2.0 Mq. 

In Figure [H we present the distributions of the physical and orbital parameters determined 
from modeling. The histograms are based on the chains obtained from MCMC running, where the 
dark and light shaded ones are with and without the constraint from blended flux, i.e. M < 2.0Mq, 
respectively. We measure the physical and orbital parameters and their uncertainties as the mean 
and standard deviation of the values in the chains and list them in Table [2j It is found that 
OGLE-2005-BLG-018 was produced by a binary lens located in the Galactic bulge with a distance 
to the lens of -Dl = 6.7 it 0.3 kpc. The lens is composed of two main-sequence stars with masses 
Ml = 0.9 lb 0.3 Mq and M2 = 0.5 it 0.1 Mq. The mass of the lens system is consistent with the 
restriction of M = M1 + M2 < 2 Mq, that was given by the blended flux. The two lens components 
are separated by a semi- major axis of a = 2.5 zt 1.0 AU and orbiting each other with an orbital 
period of P = 3.1 ± 1.3 yr. 



5. Conclusion 



We analyzed the light curve of OGLE-2005-BLG-018 based on the combined data from survey 
and follow-up observations. The light curve shows noticeable deviations from the best-fit model 
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Fig. 5. — Geometry of the binary-lens system responsible for the lensing event OGLE-2005-BLG- 
018. The filled dots represent the locations of the lens components at the times of 3 different major 
perturbations, where the bigger dots represent the heavier lens component. The insets show the 
zoom of the lens positions. The closed figure composed of concave curves represent the caustics 
where the colors correspond to those of the lens. The line with an arrow represents the source 
trajectory. The coordinates (^,r/) are centered at the center of mass of the binary and all lengths 
are scaled by the Einstein radius corresponding to the total mass of the binary lens. 

based on the standard binary parameters. Prom modeling including various higher-order effects, 
we found that the dominant second-order effect for the deviation is the orbital motion of the lens. 
Based on modeling considering full-Keplerian orbital motion and the parallax effect, we were able to 
measure the physical and orbital parameters of the lens system. Detections of higher-order effects 
and determinations of the physical lens parameters were possible due to the well-resolved light curve 
covering all three major perturbations. Unfortunately, events with such detections of higher-order 
effects are rare for events detected in current lensing experiments based on the survey /follow- up 
mode. This is because it is difficult to resolve perturbations from survey observations alone and 
even if perturbations are detected and an alert is issued for follow-up observations, it is unavoidable 
to miss part of the perturbation due to the time gap between the alert and the initiation of follow-up 
observations. 

However, the advent of next-generation experiments based on ground and in space will make it 
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possible to routinely measure higher-order effects for a large fraction of lensing events. The OGLE 
and MOA survey groups have recently upgraded their cameras with wider field of view. The Korea 
Microlensing Telescope Network (KMTNet) is a funded project that plans to achieve ~ 10 minute 
sampling of all lensing events using a network of 1.6 m telescopes to be located in three different 
continents of southern hemisphere with cameras of 4 deg^ field of view. Furthermore, t here are 
plann ed lensing surveys in space including EUCLID (iBeaulieu et al.ll201Cll ) and WFIRST (iBennett 



2010 ). that are proposed to ESA and recommended as the top ranked large space missions of 



NASA for the next decade, respectively. With these experiments come on line, nearly all events 
will be densely observed, making it possible to routinely measure the higher-order effects and thus 
to constrain the physical parameters of lenses. 



This work was supported by Creative Research Initiative program (2009-0081561) of National 
research Foundation of Korea. Work by A.G. was supported by NSF grant AST-0757888. The 
OGLE project has received funding from the European Research Council under the European 
Community's Seventh Framework Programme (FP7/2007-2013) / ERG grant agreement no. 246678 
to AU. 
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Fig. 6. — Histograms of the physical and orbital parameters, The light and dark shaded histograms 
are with and without the constraint from blended flux, respectively. 



Table 1: Best-fit Parameters for 6 Tested Models 



parameters 


standard 


parallax 


xallarap 


model 

orbit 
(2-parameters) 


parallax + 
orbit (2 parameters) 


parallax + 
orbit (4 parameters) 


Y^/dof 


3465.1/(1602) 


2352.2/(1600) 


2314.0/(1597) 


1611.8/(1600) 


1610.506/(1598) 


1607.121/(1596) 


to (HJD-2450000) 


3514.565±0.007 


3514.577±0.008 


3514.546±0.011 


3514.931±0.014 


3514.906zb0.0143 


3514.927±0.016 




0.124±0.001 


0.122±0.001 


0.126±0.001 


0.127±0.001 


0.127±0.001 


0.128±0.001 


tp, (days) 


50.44±0.09 


51.629±0.072 


49.921±0.173 


52.297±0.083 


52.337±0.115 


52.130±0.159 


S 1 


0.715±0.001 


0.702±0.001 


0.705±0.001 


0.722±0.001 


0.722±0.001 


0.724±0.001 





0.521±0.001 


0.528±0.002 


0.555±0.004 


0.539±0.003 


0.536±0.003 


0.539±0.002 


a (rad) 


4.998±0.001 


5.002±0.001 


4.998±0.002 


5.028±0.001 


5.025±0.002 


5.026±0.002 




0.025±0.001 


0.025±0.001 


0.025±0.001 


0.025±0.001 


0.025±0.001 


0.026±0.001 


" Izj,iV 




0.115±0.011 






-0.044±0.030 


-0.011±0.028 






0.342±0.008 






-0.006±0.010 


0.021±0.010 








-0.039±0.004 














-0.039±0.001 














3.98±0.25 














1.50±0.10 








J^s (yr) 






0.45±0.01 
















-0.409±0.009 


-0.387±0.011 


-0.389±0.013 


da/dt (yr-i) 








0.272±0.022 


0.328±0.049 


0.315±0.046 














-0.832±0.180 


(yr-1) 












0.581±0.161 



Notes. HJD'=HJD-2450000. 
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Table 2. Physical and Orbital Parameters 



parameter 


values 


Mtotal (Mq) 


1.38±0.39 


Ml (Mq) 


0.90±0.25 


M2 (Mo) 


0.48±0.14 


Dl (kpc) 


6.74±0.32 


a (AU) 


2.46±0.97 


P (yr) 


S.lOil.SO 


e 


0.97±0.01 


i (deg) 


-55.01±6.69 


tperi (HJD') 


2670ib352 



